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Abstract 

In this paper, a game-theoretic model for studying power control for wireless data networks in 
frequency-selective multipath environments is analyzed. The uplink of an impulse-radio ultrawideband 
system is considered. The effects of self-interference and multiple-access interference on the performance 
of generic Rake receivers are investigated for synchronous systems. Focusing on energy efficiency, a 
noncooperative game is proposed in which users in the network are allowed to choose their transmit 
powers to maximize their own utilities, and the Nash equilibrium for the proposed game is derived. 
It is shown that, due to the frequency selective multipath, the noncooperative solution is achieved at 
different signal-to-interference-plus-noise ratios, depending on the channel realization and the type of 
Rake receiver employed. A large-system analysis is performed to derive explicit expressions for the 
achieved utilities. The Pareto-optimal (cooperative) solution is also discussed and compared with the 
noncooperative approach. 
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I. Introduction 

As the demand for wireless services increases, the need for efficient resource allocation and interference 
mitigation in wireless data networks becomes more and more crucial. A fundamental goal of radio resource 
management is transmitter power control, which aims to allow each user to achieve the required quality 
of service (QoS) at the uplink receiver without causing unnecessary interference to other users in the 
system. Another key issue in wireless system design is energy consumption at user terminals, since, in 
many scenarios, the terminals are battery-powered. 

Recently, game theory has been used as an effective tool to study noncooperative power control in data 
networks [l]-[8]. The advantages of noncooperative (distributed) approaches with respect to a cooperative 
(centralized) approach are mainly due to the scalability of the network. Many of the problems to be solved 
in a communications system are in fact known to be NP-hard. As a consequence, real-time solution of 
these optimization problems in a centralized fashion becomes infeasible as the network size increases 
and as the number of users varies [1]. Game theory is the natural framework for modeling and studying 
these kinds of interactions. A model for power control as a noncooperative game is proposed in [2]. In 
this scenario, the users choose their transmit powers to maximize their utilities, defined as the ratio of 
throughput to transmit power. In [3], a network-assisted power-control scheme is proposed to improve 
the overall utility of a direct-sequence code-division multiple access (DS-CDMA) system. In [4], [5], the 
authors use pricing to obtain a more efficient solution for the power control game. Joint network-centric 
and user-centric power control are discussed in [6]. In [7], the authors propose a power control game for 
multicarrier CDMA (MC-CDMA) systems, while in [8] the effects of the receiver have been considered, 
particularly extending the study to multiuser detectors and multiantenna receivers. 

This work considers power control in ultrawideband (UWB) systems. UWB technology is considered 
to be a potential candidate for next-generation short-range high-speed data transmission, due to its 
large spreading factor (which implies large multiuser capacity) and low power spectral density (which 
allows coexistence with incumbent systems in the same frequency bands). Commonly, impulse-radio (IR) 
systems, which transmit very short pulses with a low duty cycle, are employed to implement UWB systems 
[9], [10]. In an IR system, a train of pulses is sent and the information is usually conveyed by the position 
or the polarity of the pulses, which correspond to pulse position modulation (PPM) and binary phase 
shift keying (BPSK), respectively. To prevent catastrophic collisions among different users and, thus, to 
provide robustness against multiple access interference (MAI), each information symbol is represented 
by a sequence of pulses; the positions of the pulses within that sequence are determined by a pseudo- 
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random time-hopping (TH) sequence that is specific to each user [10]. In "classical" impulse radio, the 
polarity of those pulses representing an information symbol is always the same, whether PPM or BPSK 
is employed [10], [11]. Recently, pulse-based polarity randomization was proposed [12], where each 
pulse has a random polarity code in addition to the modulation scheme, providing additional robustness 
against MAI [13] and helping to optimize the spectral shape according to US Federal Communications 
Commission (FCC) specifications [14]. Due to the large bandwidth, UWB signals have a much higher 
temporal resolution than conventional narrowband or wideband signals. Hence, channel fading cannot 
be assumed to be flat [15], and self-interference (SI) must be taken into account [16]. To the best of 
our knowledge, this paper is the first to study the problem of radio resource allocation in a frequency- 
selective multipath environment using a game-theoretic approach. Previous work in this area has assumed 
flat fading [17]-[19]. 

Our focus throughout this work is on energy efficiency. In this kind of applications it is often more 
important to maximize the number of bits transmitted per Joule of energy consumed than to maximize 
throughput. We first propose a noncooperative (distributed) game in which users are allowed to choose 
their transmit powers according to a utility-maximization criterion. Focusing on Rake receivers [20] 
at the base station, we derive the Nash equilibrium for the proposed game, also proving its existence 
and uniqueness. Using a large-system analysis, we obtain an approximation of the interference which is 
suitable for any kind of fading channels, including both small- and large-scale statistics. We also compute 
explicit expressions for the utility achieved at the equilibrium in a particular scenario and compare the 
performance of our noncooperative approach with the optimal cooperative (centralized) solution. It is 
shown that the difference between these two approaches is not significant for typical networks. 

The remainder of this paper is organized as follows. Some background for this work is given in 
Sect. JI] In Sect. [Till we describe our power control game and analyze the Nash equilibrium for this 
game. In Sect. |IVJ we use the game-theoretic framework along with a large-system analysis to evaluate 
the performance of the system in terms of transmit powers and achieved utilities. The Pareto-optimal 
(cooperative) solution to the power control game is discussed in Sect. |VJ and its performance is compared 
with that of the noncooperative approach. Numerical results are discussed in Sect. |VIJ where also an 
iterative and distributed algorithm for reaching the Nash equilibrium is presented. Some conclusions are 
drawn in Sect. I VII I 
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II. Problem Statement 

A. Formulation 

Consider the uplink of an IR-UWB data network, where every user wishes to locally and selfishly 
choose its action to maximize its own utility function. The strategy chosen by a user affects the perfor- 
mance of the other users in the network through MAI. Furthermore, since a realistic IR-UWB transmission 
takes place in frequency-selective multipath channels, the effect of SI cannot be neglected. 

To pose the power control problem as a noncooperative game, a suitable definition of a utility function 
is needed to measure energy efficiency for wireless data applications. A tradeoff relationship exists 
between obtaining high signal-to-interference-plus-noise ratio (SINR) levels and consuming low energy. 
These issues can be quantified [2] by defining the utility function of the kth user to be the ratio of its 
throughput Tfc to its transmit power p^, i.e. 

"Jfe(p) = — , (!) 
Pk 

where p = [pi, . . . ,Pk] is the vector of transmit powers, with K denoting the number of users in the 
network. Throughput, here referred to as the net number of information bits that are received without 
error per unit time, can be expressed as 

T k = j^Rkfsilk), (2) 

where D and M are the number of information bits and the total number of bits in a packet, respectively; 
Rk and 7^ are the transmission rate and the SINR for the kth user, respectively; and / s (7fc) is the efficiency 
function representing the packet success rate (PSR), i.e., the probability that a packet is received without 
an error. Our assumption is that a packet will be retransmitted if it has one or more bit errors. The 
PSR depends on the details of the data transmission, including its modulation, coding, and packet size. 
To prevent the mathematical anomalies described in [2], we replace PSR with an efficiency function 
f(lk) when calculating the throughput for our utility function. In the case of BPSK TH-IR systems in 
multipath fading channels, a reasonable approximation to the PSR for moderate-to-large values of M is 
/(Tfc) = (1 - e-W2)Af. 

However, our analysis throughout this paper is valid for any efficiency function that is increasing, S- 
shapedQand continuously differentiable, with /(0) = 0, /(+00) = 1, and /'(0) = df (jk) / djk\- /k =o = 0. 

'An increasing function is S-shaped if there is a point above which the function is concave, and below which the function is 
convex. 
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Fig. 1. User's utility as a function of transmit power for fixed interference. 



These assumptions are valid in many practical systems. Furthermore, we assume that all users have the 
same efficiency function. Generalization to the case where the efficiency function is dependent on k is 
straightforward. Note that the throughput in ((2]) could also be replaced with the Shannon capacity 
formula if u^.(p) in £T|) is appropriately modified to ensure itfc(p) = when p^ = 0. 
Combining CO) and (O, and replacing the PSR with /(7fc), 

ttfe(p) = TT^fc • (3) 

This utility function, which has units of bits/Joule, represents the total number of data bits that are 
delivered to the destination without an error per Joule of energy consumed, capturing the tradeoff between 
throughput and battery life. For the sake of simplicity, we assume that the transmission rate is the same 
for all users, i.e., R\ = ■ ■ ■ = Rk = R- All the results obtained here can easily be generalized to the 
case of unequal rates. Fig. Q] shows the shape of the utility function in © as a function of transmit power 
keeping other users' transmit power fixed (the meaning of p* and u* will be provided in the remainder 
of the paper). 
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B. System Model 

We consider a BPSK random TH-IR system^ with K users in the network transmitting to a receiver 
at a common concentration point. The processing gain of the system is assumed to be N = Nf ■ N c , 
where Nf is the number of pulses that represent one information symbol, and N c denotes the number 
of possible pulse positions in a frame [9]. The transmission is assumed to be over frequency selective 
channels, with the channel for user k modeled as a tapped delay line: 

L 

c k (t) = «i (fc) <*(< - ( l - l ) T c -n), (4) 
1=1 

where T c is the duration of the transmitted UWB pulse, which is the minimum resolvable path interval; 

(k) (k)tT 

L is the number of channel paths; otk = [a\ , a L J and are the fading coefficients and the 
delay of user k, respectively. Considering a chip-synchronous scenario, the symbols are misaligned by 
an integer multiple of the chip interval T c : = AfcT c , for every k, where is uniformly distributed 
in {0, 1, . . . ,N — 1}. In addition we assume that the channel characteristics remain unchanged over a 
number of symbol intervals [16]. 

Especially in indoor environments, multipath channels can have hundreds of multipath components 
due to the high resolution of UWB signals. In such cases, linear receivers such as matched filters (MFs), 
pulse-discarding receivers [21], and multiuser detectors (MUDs) [22] cannot provide good performance, 
since more collisions will occur through multipath components. In order to mitigate the effect of multipath 
fading as much as possible, we consider a base station where K Rake receivers [20] are used[] The Rake 
receiver for user k is composed of L fingers, where the vector (3 k = G-a^ = [p[ k \ . . . , fl^] T represents 
the combining weights for user k, and the LxL matrix G depends on the type of Rake receiver employed. 
In particular, if G = I, an all-Rake (ARake) receiver is considered. 

The SINR of the kth user at the output of the Rake receiver can be well approximated^ by [16] 

h k Pk ^ 
Ik = f , (5) 



3=1 
3+k 

where a 2 is the variance of the additive white Gaussian noise (AWGN) at the receiver, and the gains are 



2, 



Throughout all the paper, we consider IR-UWB systems with polarity code randomization [12]. 
3 For ease of calculation, perfect channel estimation is considered throughout the paper. 
4 This approximation is valid for large Nf (typically, at least 5). 
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have been introduced for convenience of notation. 

By considering frequency selective channels, the transmit power of the kth user, p^, does appear not 
only in the numerator of d5), but also in the denominator, owing to the SI due to multiple paths. In the 
following sections, we extend the approach of game theory to multipath channels, accounting for the SI 
in addition to MAI and AWGN. The problem is more challenging than with a single path since every 
user achieves a different SINR at the output of its Rake receiver. 



III. The Noncooperative Power Control Game 

In this section, we propose a noncooperative power control game (NPCG) in which every user seeks 
to maximize its own utility by choosing its transmit power. Let G = [IC, {Pk}, {^fc(p)}] t> e the proposed 



noncooperative game where K, = {1, . . . ,K} is the index set for the terminal users; Pk 



P k ,Pk 



is 



the strategy set, with p and p k denoting minimum and maximum power constraints, respectively; and 
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Uk(p) is the payoff function for user k [5]. Throughout this paper, we assume p = and p fe = p > 
for all fcE/C. 

Formally, the NPCG can be expressed as 

max Ufc(p) = max u k (pk,P-k) for k = 1, . . . , K , (12) 

where p_& denotes the vector of transmit powers of all terminals except terminal k. The latter notation is 
used to emphasize that the kth user has control over its own power p k only. Assuming equal transmission 
rate for all users, (fT2l can be rewritten as 

f(lk(Pk,P-k)) f , , rs 
max for k = 1, . . . ,K, (13) 

p k eP k pk 

where we have explicitly shown that j k is a function of p, as expressed in Q. 

The solution that is most widely used for noncooperative game theoretic problems is the Nash equi- 
librium [1]. A Nash equilibrium is a set of strategies such that no user can unilaterally improve its own 
utility. Formally, a power vector p = \p\, . . . ,px] is a Nash equilibrium of G = [/C, {Pk}, {u k (p)}] if, 
for every kefC, u k (p k ,P-k) > Uk(p' k ,P-k) for all p' k G P k . 

The Nash equilibrium concept offers a predictable, stable outcome of a game where multiple agents 
with conflicting interests compete through self-optimization and reach a point where no player wishes 
to deviate. However, such a point does not necessary exist. First, we investigate the existence of an 
equilibrium in the NPCG. 

Theorem 1: A Nash equilibrium exists in the NPCG G = [JC, {P k }, {u k {p)}]- Furthermore, the 
unconstrained maximization of the utility function occurs when each user k achieves an SINR 7^ that is 
a solution of 

f'(i*kHa-i*kho,k) = fh*k), (14) 

where 

/, (SP) (R H \2 

7o,fc = -tftT = N ' 7^ — z nn - 1 ( 15 ) 

4 S) ||*-(Bf -a* + Af -P k )\\ 2 

and fill) = dfhk)/d-y k \ 7k=rk . 

The proof of Theorem [T] can be found in App. II-AI 

The Nash equilibrium can be seen from another point of view. The power level chosen by a rational 
self-optimizing user constitutes a best response to the powers chosen by other players. Formally, terminal 
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k's best response r& : P-k — > Pk is the correspondence that assigns to each G P-k the set 

rjfc(p-k) = G Pfe : u k (p k ,P- k ) > u k (p k ,v-k)} 

for all p' fe G P k , (16) 

where P_£ is the strategy space of all users excluding user k. 

With the notion of a terminal's best response correspondence, the Nash equilibrium can be restated in 
a compact form: the power vector p is a Nash equilibrium of the NPCG G = [fC, {Pk}, { u k(p)}] if and 
only if pk G rfc(p_&) for all k G K. 

Prop. 1: Using the above definition in the NPCG, with a slight abuse of notation, terminal fc's best 
response to a given interference vector p_^ is [5] 

rjfe(p-jfe) = min(p,p£) 3 (17) 

where 

p% = arg maxu k (pk,P-k) 



(18) 



4 SP) (i - r k M 

is the unconstrained maximizer of the utility in © (see Fig. [TJ. Furthermore, p* k is unique. The proof of 
Prop. Q] can be found in App. II-Bi 

It is worth noting that, at any equilibrium of the NPCG, a terminal either attains the utility maximizing 
SINR 7^ or it fails to do so and transmits at maximum power p. 

Theorem 2: The NPCG has a unique Nash equilibrium. 
The proof of Theorem [2] can be found in App. II-CI 

IV. Analysis of the Nash equilibrium 

In the previous section, we have proven that a Nash equilibrium for the NPCG exists and is unique. 
In the following, we study the properties of this equilibrium. It is worth emphasizing that, unlike the 
previous work in this area, 7^ is dependent on k, because of the SI in ([5]). Hence, each user attains 
a different level of SINR. More importantly, the only term dependent on k in (fl4l) is 70^, which is 
affected only by the channel of user k. This means that 7? can be assumed constant when the channel 
characteristics remain unchanged, irrespectively of the transmit powers p and the channel coefficients of 
the other users. For convenience of notation, we can express 7? as a function of 70^: 

il = r (70,*) • (19) 
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2 4 6 8 10 12 14 16 18 20 



% [dB] 

Fig. 2. Shape of 7^ as a function of 70^ (M = 100). 



Fig. [2] shows the shape of 7^ as a function of 70^, where the efficiency function is taken as f{^k) = 
(1 — e _7fc//2 ) M , with M = 100. Even though 7^ is shown for values of 70^ approaching OdB, it is worth 
emphasizing that 70^ > 10 dB in most practical situations. 

As can be noticed, the NPCG proposed herein represents a generalization of the power control games 
discussed thoroughly in literature [2]-[8]. If L = 1, i.e. in a flat-fading scenario, we obtain from (|7]) and 
( fT51 ) that 7 0> fe = 00 for all k. This implies that 7^ is the same for every k € K,, and thus it is possible to 
apply the approach proposed, e.g., in [5]. 

Assumption 1: To simplify the analysis, let us assume the typical case of multiuser UWB systems, 
where N S> K. In addition, p is considered sufficiently large that Pk < P for those users who achieve 
7£. In particular, when N S> K, at the Nash equilibrium the following property holds: 

4 SP Vfc ^ Q > Vfc G £. (20) 

The heuristic derivation of (l20l) can be justified by SI reduction due to the hypothesis N ^> K > 1. 
Using CO), 7o,fc S> 1 for all k. Hence, the noncooperative solution will be similar to that studied, e.g., in 
[8]. The validity of this assumption will be shown in Sect. IVII through simulations. 

The following proposition helps identify the Nash equilibrium for a given set of channel realizations. 
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Prop. 2: A necessary and sufficient condition for a desired SINR 7^ to be achievable is 

it'fal + G 1 ) <i vfce/c, (2i) 

where 70 k is defined as in ([TBI ), and 



-K" l(MAI) 



c x = Ew- (22) 

j= i h- 

When (1211 ) holds, each user can reach the optimum SINR, and the minimum power solution to do so is 
to assign each user k a transmit power 

„• - _J ^"t ,23) 

*"*T 1 - t; -(^ + <,:')' 

When (1211) does not hold, the users cannot achieve 7^ simultaneously, and some of them would end up 
transmitting at the maximum power p. The proof of Prop. [2] can be found in App. II-DI 

Based on Prop. [2j the amount of transmit power pt required to achieve the target SINR 7^ will depend 
not only on the gain h^ F \ but also on the SI term (through 70^) and the interferers h^*^ (through 
Cfc). In order to derive some quantitative results for the utility function and for the transmit powers 
independent of SI and MAI terms, it is possible to resort to a large systems analysis. For convenience of 
notation, we introduce the following definitions, with Var[-] denoting the variance of a random variable: 

• let D" be a diagonal matrix whose elements are 



{D»}; = VVar[af]; (24) 

• let be a diagonal matrix whose elements are 

{Df } ; = yVar [/f ]; (25) 

• let C° be an L x (L — 1) matrix whose elements are 



{c?} „ = (26) 



• let be an L x (L — 1) matrix whose elements are 



let ip(-) be the matrix operator 
where Tr(-) is the trace operator. 



{C&^/^M; (27) 



¥>(•)= lim y ■&(•), (28) 

L^OO Li 
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Theorem 3: Assume that a[ are zero-mean random variables independent across k and I, and G is 
a deterministic diagonal matrix (thus implying that a k and fy- are dependent only when j = k and 
m = I). In the asymptotic case where K and Nf are finite^ while L,N C — > oo, the term CT 1 converges 
almost surely (a.s.) to 



Theorem 4: Assume a? and G as in Theorem [3] In the asymptotic case where K and iVj are finite, 
while L, N c — > oo, 

1 L-l i 

-1 a.s. 1 i=l 1=1 , an , 

7 °' fc TV n — ^v? ' m 

H D £ D *)) 

where ^ is defined as in (fTTT) and 

fc (i,L + Z-*) = {Dg},{D2}L+i-< 

+ {Df},{Dg} i+ ,_ i . (31) 
The proofs of Theorems [3] and |4] can be found in App. ITT] 

It is worth emphasize that the results above can be applied to any kind of fading models, since only 
the second-order statistics are required. Furthermore, due to the symmetry of d29l and (l30l) . it is easy to 
verify that the results are independent of large-scale fading models. Hence, Theorems [3] and |4] apply to 
any kind of channel, which may include both large- and small-scale statistics. 

For ease of calculation, in the following we derive the asymptotic values when considering a flat 
Power Delay Profile (PDP) [23] for the channel coefficients. In addition, the variance of a[ is assumed 
dependent on the user k, but independent of I, i.e., Var[a^ ] = a\ for all I. 

Prop. 3: Under the above mentioned hypotheses, when using an ARake, and thus G = I, 

- (32) 



1 a^s. V (p) 

%,k - — (33) 

5 In order for the analysis to be consistent, and also considering regulations by the FCC [14], it is worth noting that Nf could 
not be smaller than a certain threshold (Nf > 5). 
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where p = N c /L, < p < oo, and 

( 2 

1/00 = 



., v 3 - 3p + p 2 ) , p < 1, 

3 (34) 



2/(3p), p>l. 

Using definitions (fTTT) and (l24l - d28l ), and applying Theorems [3] and [4[ after some algebraic manipulations, 
the proof is straightforward. As already noticed to justify Ass. [T] but also from d33l ), "fo,k 3> 1 f° r an 
Thus, 7^ = T (70,/c) approaches 7* = T (00), leading to a nearly SINR-balancing scenario. 

An immediate consequence of Prop. [3] is the expression for transmit powers and utility functions at 
the Nash equilibrium, which are independent of the channel realizations of the other users, and of the 

si@ 



< * hr ■ TjRk^^ ' y " ni J , (36) 



,(sp) D_ f( r )(l-j*.[K-l + v(p)]/N) 
M a 2 T 

where ut is the utility function of the user k at the Nash equilibrium (see Fig. Q]), and where the condition 
(I2T1) translates into 

K-l + v(pY 



N f > 



7 



(37) 



N c 

The validity of the claims above is verified in Sect. [VI] using simulations. To show that the results (|32l ) 
- (l33l) represent good approximations not only for the model with flat PDP and equal variances, which 
has been used only for convenience of calculation, simulations are carried out using the exponential 
decaying PDP described in [24], which provides a more realistic channel model for the UWB scenario. 

V. Social Optimum 

The solution to the power control game is said to be Pareto-optimal if there exists no other power 
allocation p for which one or more users can improve their utilities without reducing the utility of any 
of the other users. It can be shown that the Nash equilibrium presented in the previous section is not 
Pareto-optimal. This means that it is possible to improve the utility of one or more users without harming 
other users. On the other hand, it can be shown that the solution to the following social problem gives 
the Pareto-optimal frontier [8]: 

K 

p opt = arg max^£ fc -u fc (p), (38) 



fc=i 



6 Of course, the amount of transmit power pi needed to achieve jl is dependent on the channel realization. 
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for £fc G E + (the set of positive real numbers). Pareto-optimal solutions are, in general, difficult to obtain. 
Here, we conjecture that the Pareto-optimal solution occurs when all users achieve the same SINRs, ^ pt- 
This approach is chosen not only because SINR balancing ensures fairness among users in terms of 
throughput and delay [8], but also because, for large systems, the Nash equilibrium is achieved when all 
SINRs are similar. We also consider the hypothesis £1 = • • • = £x = 1, suitable for a scenario without 
priority classes. Hence, the maximization d38l ) can be written as 



p opt = arg max / (7) V) — . (39) 

In a network where the hypotheses of Ass. [H Theorem [3] and Theorem 0] are fulfilled, and where 
ARake receivers are employed, at the Nash equilibrium all users achieve a certain output SINR 7 with 



h\ 'pk ^9(7). where 
with p = N c /L. Therefore, d39l can be expressed as 



<t 2 7 

1-7- A -1 + 1/ {p)\/N 



7o P t - arg max — — } h\ (41) 

since there exists a one-to-one conespondence between 7 and p. It should be noted that, while the 
maximizations in (TT3l consider no cooperation among users, (f39b assumes that users cooperate in choosing 
their transmit powers. That means that the relationship between the user's SINR and transmit power will 
be different from that in the noncooperative case. 

Prop. 4: In a network where L, N c — ► 00 and N 3> K, using ARake receivers, the Nash equilibrium 
approaches the Pareto-optimal solution. 

Proof: The solution 7^ to (|4TT ) must satisfy the condition d (/(7)/g (7)) /d7l7=7 opt = 0- Using this 
fact, combined with d40l ), gives us the equation that must be satisfied by the solution of the maximization 
problem in (HTI ): 



f'iiopthopt 



K — l + v(p) 

1 - lopt ■ - 



f inapt)- (42) 



N 

We see from (|42b that the Pareto-optimal solution differs from the solution (O of the noncooperative 
utility-maximizing method, since ((42)) also takes into account the contribution of the interferers. In 
particular, 
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Since the function T (•) is increasing with its argument for any S-shaped 7(7) (as can also be seen in 
Fig. |U), and since N/[(K-l + v (p))] < 7 , fc for all k (from © and (O), 

lopt <7<7*, (44) 

due to (fT9l ) and (l43l) . On the other hand, assuming N ^> K and < p < 00 implies "y opt —> 7*. From 
(l44l ). it is apparent that 7 — > 7* as well. This means that, in almost all typical scenarios, the target 
SINR for the noncooperative game, 7, is close to the target SINR for the Pareto-optimal solution, ") pt- 
Consequently, the average utility provided by the Nash equilibrium is close to the one achieved according 
to the Pareto-optimal solution. ■ 
The validity of the above claims will be verified in Sect. |VI] using simulations. 

VI. Simulation Results 

A. Implementation 

In this subsection, we present an iterative and distributed algorithm for reaching a Nash equilibrium 
of the proposed power control game. This algorithm is applicable to all types of Rake receivers, as well 
as to any kind of channel model. The description of the algorithm is as follows. 



The Best-Response Power-Control (BRPC) Algorithm 



Consider a network with K users, a processing gain N = Nf ■ N c , a channel with L fading paths, 
and a maximum transmit power p. 

1) Simulate the channel fading coefficients atk for all users according to the chosen channel model. 

2) Set the Rake receivers coefficients (3 k for all users according to the chosen receiver. 

3) Compute the SP term, hf P \ the SI term, hf l) , the MAI term, h { ™ M) , according to © - ®, and 
the optimum SINR, 7^, solution of ([141) . for all users. 

4) Initialize randomly the transmit powers of all users within the range [0,p]. 

5) Compute the received SINR 7^ at the base station for each user according to ([5]). 

6) Set k = 1. 

7) Adjust the transmit power for user k according to (fT71) and to 

Pk-Pk ~, —/ , (43) 



where 70 k is defined as in ( |T5l ). 
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8) k = Jfe + l. 

9) If k < K, then go back to Step 7. 

10) Stop if the powers have converged; otherwise, go to Step 5. 

This is a best-response algorithm, since at each stage a user decides to transmit at a power that 
maximizes its own utility (i.e., its best-response strategy), given the current conditions of the system. 
Looking at Step 7, it may appear from ( fT8l ) that the each user should know its own transmit power p^ 
and ratio 70^, as well as some other quantities (jjj and h'^ A1 for j ^ k) relevant to all of the other users 
in the network. On the contrary, it turns out that user k only needs to know its own received SINR at 
the base station 7^. In fact, the term due to interference-plus-noise in ( fT8l ) can be obtained from (f5]) as 



where pt is the transmit power of user k at the previous iteration. Therefore, after straightforward 
manipulation, (TT8T ) translates into the noncooperative update (l45T ). The received SINR 7fc can be fed back 
to the user terminal from the access point, along with SP and SI terms. 

It is clear that the above algorithm converges to a Nash equilibrium, whose existence and uniqueness 
are proven in App. U The convergence of the BRPC algorithm is also validated in Sect. IVI-BI using 
extensive simulations. 

B. Numerical Results 

In this subsection, we present numerical results for the analysis presented in the previous sections. 
We assume that each packet contains 100b of information and no overhead (i.e., D = M = 100). The 
transmission rate is R = lOOkb/s, the thermal noise power is a 2 = 5 x 10~ 16 W, and the maximum 
transmit power is p = 1 iiW. We use the efficiency function /(7&) = (1 — e _7fc / 2 ) M . Using M = 100, 
7* = r (00) = 11.1 dB. To model the UWB scenario [15], channel gains are simulated following [24], 
where both small- and large-scale statistics are taken into account. The distance between the users and 
the base station is assumed to be uniformly distributed between 3 and 20 m. 

Before showing the numerical results for both the noncooperative and the cooperative approaches, 
some simulations are provided to verify the validity of Ass. [j] introduced in Sect. JV] Table U reports 
the ratio cr 2 /rj 2 of the variance a 2 to the squared mean value rj 2 of the values q = h^ V ^p* k , obtained 
averaging 10 000 realizations of channel coefficients for different values of network parameters using 
ARake receivers. We can see that, when the processing gain is much greater than the number of users, 
(Tq/r] 2 <C 1. Hence, (1201 can be used to carry out the theoretical analysis of the Nash equilibrium. 
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TABLE I 

Ratio o^/t^ for different network parameters. 





(L,K) 


(N c ,N f ) 


(20,8) 


(20,16) 


(50,8) 


(50,16) 


(30,10) 


9.4E-4 


3.2E-3 


4.8E-4 


1.7E-3 


(30,50) 


2.9E-5 


6.4E-5 


1.6E-5 


3.4E-5 


(50,10) 


2.9E-4 


6.8E-4 


1.5E-4 


3.7E-4 


(50,50) 


1.0E-5 


2.2E-5 


0.6E-5 


1.2E-5 


(100,10) 


6.7E-5 


1.5E-4 


3.7E-5 


7.8E-5 


(100,50) 


0.3E-5 


0.6E-5 


0.1E-5 


0.3E-5 



Fig. [3] shows the utilities achieved at the Nash equilibrium as functions of the channel gains = 
||afc|| . These results have been obtained using random channel realizations for K = 16 users. The 
number of possible pulse positions is 7Y C = 100, while the number of paths is L = 60, in order to 
satisfy the large system assumption with v (p) = 0.4. The number of frame is Nf = 10, thus leading 
to a processing gain N = 1000 K. The line represents the theoretical values of (l36l ) when using an 
ARake, whereas the square markers correspond to the simulations achieved with the BRPC algorithm. 
We can see that the simulations match closely with the theoretical results. 

Fig. |4] shows the utility as a function of the channel gain when the processing gain N is constant, but 
the ratio N c /Nf is variable^ The results have been obtained for a network with K = 16 users, L = 100 
channel paths and processing gain N = 1000, using ARake receivers at the base station. The solid 
line corresponds to N c /Nj = 40, the dotted line represents N c /Nf = 2.5, and the dashed line shows 
N c /Nf = 0.1. As expected, higher N c /Nf ratios (and thus higher iV c , when N is fixed) correspond to 
higher utility, as described in d36l) . In fact, v (p) decreases as p increases, i.e., as N c increases (when 
L is fixed). This result complies with theoretical analysis of UWB systems [16], since, for a fixed total 
processing gain N, increasing the number of chips per frame, N c , will decrease the effects of SI, while 
the dependency of the expressions on the MAI remains unchanged. Hence, a system with a higher A^ c 
achieves better performance. 

Fig. [5] shows the probability P a = Prjmaxfc p^. =p = l pW} of having at least one user transmitting 

7 To avoid a too busy graph, only the theoretical values are reported. 
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Fig. 3. Achieved utility versus channel gain at the Nash equilibrium for the ARake receiver. 



at the maximum power, as a function of the number of frames Nf. We consider 10 000 realizations of 
the channel gains, using a network with ARake receivers at the base station, K = 32 users, iV c = 50, 
and L = 100 (thus p = 0.5 and v (p) ~ 1.17). From (l37l ). the minimum value of Nf that allows all K 
users to achieve the optimum SINRs is Nf = [8.33] = 9. The simulations thus agree with the analytical 
results of Sect. |IV] 

We now analyze the performance of the system when using a Pareto-optimal solution instead of the 
Nash equilibrium. Fig. [6] shows the normalized utility u^/h^ as a function of the load factor p. We 
consider a network with K = 5 users, Nf = 20 frames and ARake receivers at the base station. 
The lines represent theoretical values of Nash equilibrium (dotted line), using d36l ), and of the social 
optimum solution (solid line), using (|36l ) again, but substituting 7* with the numerical solution of ((42]) . 
7o p4 . The markers correspond to the simulation results. In particular, the circles represent the averaged 
solution of the BRPC algorithm, while the square markers show averaged numerical results (through a 
complete search) of the maximization d38l ), with ^ = 1. As stated in Sect. |Vl the difference between 
the noncooperative approach and the Pareto-optimal solution is not significant. Fig. [7] compares the 
target SINRs of the noncooperative solutions with the target SINRs of the Pareto-optimal solutions. As 
before, the lines correspond to the theoretical values (dashed line for the noncooperative solution, solid 
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Fig. 4. Achieved utility versus channel gain at the Nash equilibrium for different ratios N c /Nf. 



line for the social optimum solution), while the markers represent the simulation results (circles for the 
noncooperative solutions, square markers for the Pareto-optimal solution). It is seen that, in both cases, 
the averaged target SINRs for the Nash equilibrium, 7, and the averaged target SINRs for the social 
optimum solution, jopt, are very close to 7*, as shown in Sect. [V] 

VII. Conclusion 

In this paper, we have used a game-theoretic approach to study power control for a wireless data 
network in frequency-selective environments, where the user terminals transmit IR-UWB signals and 
the common concentration point employs Rake receivers. A noncooperative game has been proposed 
in which users are allowed to choose their transmit powers according to a utility-maximizing criterion, 
where the utility function has been defined as the ratio of the overall throughput to the transmit power. 
For this utility function, we have shown that there exists a unique Nash equilibrium for the proposed 
game, but, due to the frequency selective multipath, this equilibrium is achieved at a different output 
SINR for each user, depending on the channel realization and the kind of Rake receiver used. Resorting 
to a large system analysis, we have obtained a general characterization for the terms due to multiple 
access interference and self-interference, suitable for any kind of channel model and for different types 
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Fig. 5. Probability of having at least one user transmitting at maximum power versus number of frames. 

of Rake receiver. Furthermore, explicit expressions for the utilities achieved at the equilibrium have been 
derived for the case of an ARake receiver. It has also been shown that, under these hypotheses, the 
noncooperative solution leads to a nearly SINR-balancing scenario. In order to evaluate the efficiency of 
the Nash equilibrium, we have studied an optimum cooperative solution for the case of ARake receivers, 
where the network seeks to maximize the sum of the users' utilities. It has been shown that the difference 
in performance between Nash and cooperative solutions is not significant for typical values of network 
parameters. 

Appendix I 

A. Proof of Theorem [7] 

The following results are obtained from [5], [25]-[27]. 

A Nash equilibrium exists in the noncooperative game G = [fC, {Pk}, {uk(p)}] if, for all k = 1, . . . ,K: 

1) Pk is a nonempty, convex, and compact subset of some Euclidean space R , and 

2) Ufc(p) is continuous in p and quasi-concave in p k . 

Each user has a strategy space that is defined by a minimum power p k and a maximum power p k , 
and all power values in between. We also assume that p k > p k . Thus, the first condition is satisfied. 
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Fig. 6. Comparison of the normalized utility versus load factor for the noncooperative and Pareto-optimal solutions. 



Moreover, since pk > 0, it is apparent from © and ® that Ufc(p) is continuous in p. To show that the 
utility function u^(p) is quasi-concave in p^ for all k in the NPCG, it is sufficient to prove that the local 
maximum of «fc(p) is at the same time a global maximum [28], [29]. 

For a differentiable function, the first-order necessary optimality condition is given by duk(p)/dpk = 0. 
Recalling © and ©, the partial derivative of itfc(p) with respect to pk is 

= ^2 {ftikH (1 - 7kho,k) - /(7k)) , (47) 
<9p fc Mpl 

where 70 & is defined as in (fT3T ) and f'(jk) = ^/ (7fc) / '^7fc- F° r the sa ke of simplicity, we do not explicitly 
show the dependence of 7^ on p^. Since p*. > for the NPCG, we examine only positive real numbers. 
Evaluating (l47l ) at pp. = 0, we get duk(pk,P-k)/dpk = 0. Therefore, p^ = is a stationary point and 
the value of utility at this point is u&(p) = 0. If we evaluate utility in the e-neighborhood of p^. = 0, 
where e is a small positive number, we notice that utility is positive, which implies utility is increasing 
at pk = 0. Hence, pp. = cannot be a local maximum. For nonzero values of the transmit power, we 
examine the values of 7^ = 7fc(p£) such that duk{pk,P-k)/dpk\p k =p* k = 0, thus satisfying the first-order 
necessary optimality condition. In other words, we evaluate 7^ such that 

7*(l-7fc/70,fc) = /(7fc)//'(7fc), (48) 
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Fig. 7. Comparison of the target SINRs versus load factor for the noncooperative and Pareto-optimal solutions. 



as shown in (fT4l . 

We observe that the left-hand side of the above equation is a concave parabola with its vertex in 7^ = 
7o,fc/2 > 0, and d(7& (1 — 7fc/7o,fc)) /dlk\^ k =o = 1- The right-hand side is an increasing function, with 
d(/(7fc)//'(7fc)) /d7kW=o = l/i' < 1 when /'(0) = 0, where i! = min{i G N : d*/(7fc)/d7*k=o / 0}. 
Furthermore, the equation is satisfied at 7^ = 0. Therefore, there is a single value 7? that satisfies ([141 ) 
for 7^ > 0. The second-order partial derivative of the utility with respect to the power reveals that this 
point is a local maximum and therefore a global maximum. Hence, the utility function of user k is 
quasi-concave in pk for all k. The same conclusion applies also if p* k > p k for some k, even though 
7^ cannot be achieved. In fact, by applying the previous considerations to (l47l) . it is easy to verify that 
Uk(pk,P-k) is strictly increasing in 7^ G [0>7fc(Pfc) = 7fc)> which in turn, from ((5]), is strictly increasing 
in pk G [0,p£). Since p^ G \p k ,p k ], which is a subset of [0,p£), Uk{pk = PkiP-k) represents both the 
local and the global maximum of the utility function. ■ 

Lemma 1: The solution 7^ of (fT4b satisfies the condition 

< 7fc < 7b,* • (49) 
Proof: As /(7^) is an increasing function of 7^, /'(7j£) > for every 7^. Since the existence 
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of the solution is ensured by Theorem Q] and 7^ and fi^Z) are both greater than zero, the condition 



B. Proof of Prop. □ 

Using Theorem [TJ for a given interference, the SINR j k corresponds to the transmit power p k as in 
(fT8l) . Since 7^ is the unique maximizer of the utility, the correspondence between the transmit power 
and the SINR must be studied. As can be verified, (fT8l ) represents the equation of a hyperbola passing 
through the origin, with the asymptotes parallel to the Cartesian axes. In particular, the vertical asymptote 
is 7^ = 7q£. Therefore, using Lemma Q] presented in App. II- Al there exists a one-to-one correspondence 
between the transmit power, p* k € [0, +00), and the SINR, 7^ € [0,70^). Thus, the transmit power pt 
is also unique. If p* k £ for some user k, since it is not a feasible point, then p* k cannot be the best 
response to a given p_&. m this case, we observe that duj { (p)/dpi t < for any 7^ < 7^, and hence 
for any pf. < pt. This implies that the utility function is increasing in that region. Since p is the largest 
power in the strategy space, it yields the highest utility among all Pk <p and thus is the best response 
to p k . ■ 

C. Proof of Theorem \2\ 

By Theorem [T] we know that there exists an equilibrium in the NPCG. Let p denote the Nash 
equilibrium in the NPCG. By definition, the Nash equilibrium must satisfy p = r(p), where r(p) = 
[ri(p), . . . , 7"if (p)]. The fixed point p = r(p) is unique if the correspondence r(p) is a standard function 
[30], i.e., if it satisfies the following properties: 

1) positivity: r(p) > 0; 

2) monotonicity: if p > p', then r(p) > r(p'); 

3) scalability: for all \i > 1, /ur(p) > r(/ip). 

It is apparent that rfc(p) = rk(p~k). Taking into account (fTTT ) and (fT8T ). the first condition translates into 
p\ > for all k G K. Using ©, ® and (j49]>, the proof is straightforward. Recalling (Q/7]) and (TT8i 
the second and the third condition are also apparent, since p_^ modifies only the numerator of ( fT8j ). 
Therefore, since r(p) is a standard function, the Nash equilibrium of the NPCG is unique. ■ 

D. Proof of PropM 

Based on Prop. [T] when all users reach the Nash equilibrium, their transmit powers are 



C 1 ~ 7fc/7o,fc) > must hold. 




(50) 
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Using Ass. [TJ in (l50l) . it is straightforward to obtain: 



^-7£-(7o^ + CrjJ = ^>0, (51) 
which implies 7^ • (7^ + C^ 1 ) < 1> proving necessity. It is also straightforward to show that, if each 



terminal A; uses transmit power p* k as in d23l , all terminals will achieve the SINR requirement, finishing 
the proof of sufficiency. Finally, consider any other joint distribution of powers and channel realizations, 
and let q' = inffc g £ j^jjf^Pfc}- Then, by exactly the same argument as was used in the proof of necessity, 

4 > P\ = <!■ (52) 

This means that assigning powers according to d23l does indeed give the minimal power solution. ■ 



Appendix II 

In order to prove Theorems [3] and 01 it is worth introducing the following results. 

Lemma 2: [31] Consider an n-dimensional vector r n = [R\, . . . , R n ] with independent and identically 
distributed (i.i.d.) standardized (complex) entries (i.e. E[i2j] = and E[[i?j| 2 ] = 1, with E[-] denoting 
expectation), and let C n be an n x n (complex) matrix independent of r n . For any integer p, 

E[|rf C n r n - Tr(C n ) f] = K p ({E[\Ri l'] Tr(C n C^))^ 2 

+ E[| J R 1 |^]Tr(C n C^) p / 2 ), (53) 

where the constant K p does not depend either on n or on C n . 

Theorem 5: Consider an n-dimensional vector x n = -^[Xi, . . . ,X n ] with i.i.d. standardized (com- 
plex) entries with finite eighth moment, and let C n be an n x n (complex) matrix independent of x n 
with uniformly bounded spectral radius for all n. Under these hypotheses, 

xfC n x n a -4- lTr(C n ). (54) 
Proof: Using Lemma [2] and Markov's inequality, 

E[|x^C n x n -lTr(Cg| 4 ] 



Pr[|x^C n x n -lTr(C n )| > e] < 



e 1 



— 1% ' 9 — 7, (55) 



where k < 00 is a constant value independent of n. Thus, 

00 

Pr[|x^ C n x n - Tr(C n ) | 4 > e] < 00. (56) 

n=l 

Using the Borel-Cantelli lemma [32], the result (l54l) is straightforward. ■ 
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Theorem 6: Suppose x n = [X\ , . . . , X n ] and y„ = [Yi , . . . , Y^J are re-dimensional independent vectors 
with i.i.d. standardized (complex) entries with finite eighth moment, and C n is an n x re matrix (complex) 
independent on x n and y n with uniformly bounded spectral radius for all n. Then, 

xfC n y„ a 4-0. (57) 
Proof: The proof can be obtained using the same steps as that of Theorem [5] ■ 

A. Proof of Theorem \3\ 

To prove that N£7 converges a.s. to non-random limits, we focus on the ratio 



N 



, (MAI) 
% 



+ 



Af -/3 fc 2 + |/3f-a i | 2 



(/3f • a fc ) • (/3f • a, 



l 

L 2 



Af-/3 fc +|j3f-a i | 



(58) 

It is sufficient to show that both numerator and denominator of (1581 ) converge a.s. to a non-random limit. 
Let 



w (fc) = Qyayl ^ 



(59) 



and 



(3 k = Ga k 



(60) 



where D^? and are defined as in (124b and (1251 ). respectively; the operator o denotes the Hadamard 
(element-wise) product; and the matrix G is dependent on the type of Rake receiver employed. Using 
d59l ) and ([60]), by Theorem [5J we obtain 



1 



J2 ll B " 



i=l Z=i+1 



L-l 



m=l 



lim yJ^Xi, 



(61) 



Z=l 
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where </?(•) is defined as in (l28l) and 

are independent random variables, with 

E[xd = ^{Df}/ + iEW>- 



L 

m=l 



m=l 

and 



- , Tuil)'/) 2 ; 



< oo. 



m 

m=l 



Using the weak version of the law of large numbers for non-i.i.d. random variables, 

^ L-l L-l 

1=1 1=1 

i=l m=l 

where C£ and are defined as in d26l) and (f2Tb . respectively. Similar arguments yield 



^||Af-/3 fe |r^>(DfC«C^Df 



Then applying Theorem [6j from ( f6Qb we obtain 

-/3 fc • aj -» 0, 

since /3 fc is independent of ay. Analogously, using Theorem [51 from (|59l ) we obtain 

1/3* • ^ v ( D £d£ 
Using (|65l ) - (l68i the result d29l ) is straightforward. 
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B. Proof of Theorem |4] 

In order to prove that N/jq^ converges a.s. to a non-random limit, it is sufficient to show that both 
the numerator and the denominator converge to non-random limits. Note that 

||»- (Bf -«» + Af -/3j|| 2 = 

i=i \i=i 1=1 j 



L-l 



L+l-i 



i=l 



(69) 

where 0^ (I, L + I — i) is defined as in d3TT ). 

Following the same steps as in the proof of Theorem [3l after some algebraic manipulation, it can be 
proven that 



& £*.£!IMMzi <70) 

i=l Z=l 

Using d70l) , in conjunction with d68l ), the result d30l is straightforward. ■ 
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